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We analyze, both analytically and numerically, the time-dependence of the return probability 
in closed systems of interacting particles. Main attention is paid to the interplay between two 
regimes, one of which is characterized by the Gaussian decay of the return probability, and another 
one is the well known regime of the exponential decay. Our analytical estimates are confirmed by 
the numerical data obtained for two models with random interaction. In view of these results, we 
also briefly discuss the dynamical model which was recently proposed for the implementation of a 
quantum computation. 
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As is known, in many-body systems the density of en- 
ergy levels increases extremely fast with an excitation 
energy. As a result, for highly excited states even a 
very weak interaction between particles can lead to a 
strong mixing between the unperturbed basis states, re- 
sulting in a complex structure of exact eigenstates. In 
this case one can speak about chaotic eigenstates since 
their components can be practically treated as pseudo- 
random ones. In the dynamics, this fact results in a re- 
laxation of the system to a steady-state distribution. The 
studies of complex atoms [1], multi-charged ions [2], nu- 
clei [3], Bose-Einstein condensates [4], and spin systems 
[5,6] have confirmed the dynamical origin of statistical 
laws in isolated systems (see, also, Refs. [7,8] and the 
review [9]). 

Recently, the theory of many-body chaos has been ex- 
tended to the models of a quantum computation. In par- 
ticular, it was argued that due to a very high density 
of energy levels, any kind of perturbation may lead to 
decoherence effects thus destroying the operability of a 
quantum computer. Therefore, it is of importance to 
search for the conditions when the role of chaos can be 
significantly reduced [10,11]. 

So far, the study of the many-body chaos has been 
mainly restricted by the investigation of statistical prop- 
erties of the energy spectra and eigenstates. On the other 
hand, in view of experimental applications, one needs to 
know what are the dynamical properties of quantum sys- 
tems with strongly interacting particles. Below, we an- 
alyze the time dependence of the widely used quantity, 
the so-called return probability which determines global 
properties of the dynamics. Our main interest is in the 
interplay between the exponential and Gaussian decrease 
of the return probability, in dependence on the strength 
of the interaction between particles. 



In what follows we consider physical systems described 
by the total Hamiltonian H which can written in the sep- 
arable form, 



H = Hn 



V. 



(1) 



Typically, such a representation is used when studying 
the infiuence of a perturbation y on a system governed 
by the unperturbed Hamiltonian Hq. In many physical 
applications this form of H is quite natural, reflecting a 
different nature of Hq and V. However, quite often the 
separation of H into the two parts is not well defined, 
as in the case of mean field approaches used to introduce 
"good" variables in which Hq has a relatively simple form 
in comparison with a "residual interaction" V. In what 
follows we discuss some of generic properties of the dy- 
namics of the model (1), by studying specific forms of the 
Hamiltonian H. 

It is naturally to represent H in the unperturbed basis 
|fc) of Hq. Then the total Hamiltonian is presented by 
the sum of the diagonal matrix Hq plus the perturba- 
tion matrix V with the matrix elements Vik = (Z|T^|fc). 
In order to consider the evolution of wave packets in the 
basis of Hq, one has to express exact eigenstates \a) of 
the total Hamiltonian H in terms of basis states |fc) of 
Hq, 



C^\k). 



(2) 



The coefficients give the expansion of an exact state 
in terms of the basis states (for a fixed), or the expansion 
of a basis state in terms of the exact states (for k fixed). 
In principle, the knowledge of the state matrix C" and 
the corresponding energy spectrum E" gives a complete 
information about the system. 

Of particular interest is the case when the initial state 
\E'(0) is a basis state |fco). Then the evolution of the 
function is described by the expression. 
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m = E ^nC^o \ko) expi-iE'^t). (3) 

Here and below we assume that h = 1. As one can see, 
the probability 

w, = \A,\' = \{k\^{t))\^ (4) 

to find the system at time t in the state is determined 
by the amplitude 

Ak = (A;|exp(-ii?t)|fc) =E|CfePexp(-i£;"<). (5) 

a 

Our main interest is in the return probability Wo{t) 
which is the probability to find the system at time t in 
the initial state \ko). One can see that the return proba- 
bility is determined by the expression, 

Wo{t) = Wko(.t) = \AkM\ (6) 

where 

^fco(*)=El'^fcol'expHiS"i)~ / Pko{E)exp{-iEt)dE. 

a ■' 

(7) 

Here we replaced the summation by integration that can 
be done if the number of large components C^^ is large. 

Indeed, these components strongly fluctuate around their 
mean values and quite often can be considered as pseudo- 
random quantities. In fact, this condition of a large num- 
ber of pseudo-random components in exact eigenstates 
can be used as the definition of chaos in quantum sys- 
tems (for details, sec Refs. [12,7,13]). In this case, the 
time dependence of VFo(t) is entirely determined by the 
Fourier transform of Pfco {E) = P{E, Ek„) where E is the 
energy of exact eigenstates and Ek„ is the energy corre- 
sponding to the imperturbed state |fco). This quantity is 
known in the literature as the strength function (SF) or 
local spectral density of states, 



PiE,E,„)^\Cl:>\^p{E). (8) 

Here p{E) is the density of states of the total Hamilto- 
nian H, and the average is performed over a number of 
states with energies close to E. 

To analyze generical properties of the return proba- 
bility, let us start with its behavior at small times. Ac- 
cording to the perturbation theory, one can easily get the 
general expression, 

Wo{t)^l-Alt^ (9) 

where is the variance of the strength function in the 
unperturbed energy space, determined as 

Al=Y.V,\. (10) 

k^ko 



Note that the above expression is universal in the sense 
that it is exact for any kind of the perturbation V. Prac- 
tically, the initial time scale for the perturbative expres- 
sion (9) to be valid is very small, and the main interest is 
in the time-dependence of Wq {t) beyond this time scale. 



III. TWO-BODY RANDOM INTERACTION 
MODEL 

In order to analyze the behavior of Wo{t) on a large 
time scale, we consider the model which describes a 
closed system of N Fermi-particles occupying M single- 
particle levels of energies eg- The total Hamiltonian can 
be represented in the form (1) where 

M M 

HQ = ^esalas; V= ^ ^^SlS2S3S4 «Sl«S2"s3aS4• 

S=l Sl,S2,S3,S4 = l 

(11) 

Here Hq stands for non-interacting particles, and the in- 
teraction V between the particles is expressed in terms of 
two-body matrix elements VsiS2S3S4- The many-particle 
basis |fc) of Hq is defined by the Slater determinant, 
1^) = '^si ■ ■ ■ ^In I*^)' where al. and are the creation- 
annihilation operators, and |0) is the ground state. As 
one can see, the interaction between particles is assumed 
to have a two-body nature, therefore, each many-particle 
matrix element Vik = (Z| V|fc) is a sum of a number of two- 
body matrix elements FsiS2S3S4 involving at most four 
single- particle states |s) (for details, see, for instance, 
Ref. [8]). 

Note that the approach we consider here, is also valid 
for quasi-particles that appear in the mean-field theo- 
ries. In this case Hq stands for the mean-field part of the 
Hamiltonian, and V describes a residual interaction. The 
total Hamiltonian H in the form (11) describes generic 
properties of such physical systems as complex atoms, 
nuclei, quantum dots, etc. The energies Cg in such appli- 
cations are, in fact, renormalized quasi-particle energies. 

In many realistic applications the interaction V be- 
tween particles (quasi-particles) is so strong and com- 
plicated that practically one can describe such an in- 
teraction by assuming that all two-body matrix elements 
are distributed randomly according to some distribution. 
Thus, the simplest version of the two-body random inter- 
action model (TBRI) is the Hamiltonian (11) in which all 
matrix elements V^siS2S3S4 are random Gaussian numbers 
with the zero mean and the variance ^^^^ It is inter- 
esting to note that the many-body matrix elements Vik 
are, however, weakly correlated, due to the fact that the 
same two-body matrix elements enter in different many- 
body matrix elements Vik- In general, these correlations 
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can be neglected, however, for specific observables they 
give rise to unexpected results (for details see Ref. [14]). 

Without the loss of generality one can assume that 
the single-particle spectrum has the constant mean level 
spacing, do = (es+i—Cs) = 1; here the brackets (...) stand 
for the average over random values of Cg. The number 
of many-body states increases very fast with an increase 
of the number of particles N and number M of single- 
particles states. For this reason even for a relatively small 
number of particles the exact eigenstates may consist of 
many unperturbed basis states, thus providing us with 
a possibility to use statistical methods. In particular, 
a novel approach has been developed in Rcfs. [7,12,15], 
that is based on the chaotic structure of eigenstates in 
a given basis of unperturbed many-particle states. This 
approach allows one to relate statistical properties of ex- 
act eigenstates in many-body representation directly to 
the properties of single-particle operators, such as the 
occupation number distribution of single-particle states. 

As was shown above, the time dependence of the return 
probability is entirely determined by the Fourier trans- 
form of the SF. In many applications the SF is known 
to have the Breit-Wigner (BW) form (or, the same, the 
Lorentzian) resulting from the application of the Ran- 
dom Matrix Theory, (see, for instance, Ref. [16]). The 
half-width Tq of the BW is given by the Fermi golden 
rule and in our case reads as 

ToiE) = T{E, ko) ^ 2w\Vk^po{E). (12) 

Here |Vfe(,fe|^ is the mean square value of many-body ma- 
trix elements (obtained by the average over k), and po{E) 
is the density of those states which are directly coupled to 
the basis state \ko) by the interaction V. Note that this 
density po{E) is much smaller than total density ptot{E) 
of all many-body states. This fact is manifested by a 
large number of zero matrix elements for any fixed line 
in the Hamiltonian matrix Hi^. 

It should be stressed that the above result for the BW 
form of the SF is based on a non-perturbative approach 
according to which a large number of many-body states 
are coupled by a relatively strong interaction. As a result, 
the decrease of the return probability has the exponential 
time dependence, 

Wo(t)^exp(-rot), (13) 

apart from a small time scale i < to on which the 
quadratic decrease (9) occurs. 

For a long time it was assumed that the exponential de- 
crease of the return probability is the only regime which is 
physically relevant to the dynamics of systems with many 
interacting particles. However, recently it was found that 
in many situations the form of the SF can be quite close 
to the Gaussian (see, e.g. [3]). This fact is related to a 
finite width Ay of the interaction in the energy space for 
isolated systems. For the TBRI model it was shown [17] 



that if Fo is much less than Ay, the form of the SF, is, 
indeed, the BW. However, in the other limit, Fq > Ay, 
of a very strong interaction, the leading dependence of 
Wo{t) is the Gaussian, 

Wo{t) ~ exp(-A|t2), (14) 

and occurs on a time scale < t < tc- After, for t > tc, 
the decay of Wo{t) is described by the exponential func- 
tion [18]. 

The transition from the BW to the Gaussian for the 
TBRI model has been analyzed in Ref. [19], although 
the analytical expression in a closed form is unknown. 
In order to evaluate Wo{t), in Ref. [18] a phenomeno- 
logical expression was suggested that depends on both 
parameters, Fo and A^;. The analytical expression for 
the variance A|; of the strength function can be found 
explicitly [7], 

,,2 

Al = ^N{N-l){M-N){M-N + 3). (15) 

Here [—vq, vq] is the range within which the two-body ma- 
trix elements are distributed randomly with a constant 
probability, therefore, ^^^^ = It is interesting to 

note that for Fcrmi-particlcs the variance A|; turns out 
to be independent of a specific basis state |A:o). 

Thus, in the case of a not very strong perturbation the 
decrease of the return probability is the exponential one, 
and with an increase of the interaction V, one should ex- 
pect a quite large time scale on which the Gaussian form 
(14) occurs. Let us check these predictions by making use 
of numerical data. In all our calculations we have used 
N = 6 Fermi-particles occupying M = 12 single-particle 
states; this results in the size 924 x 924 of the Hamil- 
tonian matrix. For simplicity, the initially excited state 
was taken at the center of the energy spectrum, ko = 462, 
where the density of many-body states is maximal, and 
the energy spectrum is symmetrical. 

For a relatively weak (however, non-perturbative) in- 
teraction, vq = 0.12, the data, indeed, demonstrate a 
clear exponential dependence, up to some time scale be- 
yond which the finite size effects of the Hamiltonian ma- 
trix are important, see Fig.l. The fit to the exponential 
dependence (13) gives Fq ~ 0.97. To compare with the 
analytical expression (12), one should note that for the 
TBRI model this expression is difficult to use directly, 
since the quantity Tq{E) is not well-defined by Eq.(12). 
The problem is that the density of directly coupled many- 
body states strongly changes in dependence on k. The 
rough estimate can be obtained as follows. First, one can 

2 2 

relate the term \ Vkok\ to A| in the way, A% = K\Vkok\ , 
where K is the number of non-zero elements in any line 
of the Hamiltonian matrix (which is independent of ko 
[7]). Second, the simplest estimate of the mean density 
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po is due to the average width Ay of the interaction in 
the energy space, po ~ K/Ay- As a result, one gets, 

ro«27r^. (16) 

Finally, it can be shown that the simplest estimate for 
the width of the interaction reads as Ay w 2do{M — N). 

This gives To w 1.03 which is a good result, taking into 
account the problems with the evaluation of the expres- 
sion for Fq. 
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FIG. 1. Time dependence of the return probability Wo{t) 
for the interaction strength vo = 0.12 for which the SF has 
the BW form. Full circles stand for numerical data, and the 
straight line (shifted for a better visualization) is the linear 
fit on the time scale where the linear slope is clearly seen. 
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FIG. 2. Time dependence of the return probability Wo{t) 
for the case of the Gaussian form of the strength function. 
FHiU circles correspond to the numerical data for vq — 0.5 and 
the solid curve is the analytical dependence (14) with A% 
determined by Eq.(15). 

Now let us consider another limit case of a very strong 
interaction, vq = 0.5, when the SF is quite close to the 
Gaussian. Numerical data reported in Fig. 2 manifest a 
long Gaussian decrease of the return probability. Note 
that the deviation from the Gaussian dependence to- 
wards the exponential one (linear slope in Fig. 2 after 



t w 0.4) starts for very small values of Wo(t). There- 
fore, practically the decrease of the return probability is 
described by the dependence (14). 

Finally, we conclude with the intermediate case when 
both dependencies, (14) and (13), are important for the 
description of the return probability. In Fig. 3 one can see 
that there arc two time scales. On the first time scale, 
t < tc ~ 0.3, the decrease of Wo{t) has the Gaussian 
form (14), and for t > tc it changes to the exponential 
dependence. It is now instructive to analyze the critical 
time tc which divides these two regimes. 



One can see that tc is the time to resolve the finiteness of 
the width of the interaction. If this width is very large, 
the exponential decrease starts on a small time scale. 
Contrary, for relatively small values of Ay the Gaussian 
decrease of Wo{t) starts from t = Q and lasts for a long 
time. According to this estimate, wc have tc » 0.5 which 
roughly corresponds to the data. Note that the critical 
time tc is, in fact, not a well defined quantity and can be 
determined up to some numerical factor of the order one. 
Remarkably, the time of the correspondence of the data 
to Eq.(14) turns out to be independent of the perturba- 
tion strength. 




FIG. 3. Return probability Wo{t) for the intermediate sit- 
uation with Vq = 0.25, when both regimes, the Gaussian and 
exponential ones, are characteristic of the dynamics. Full cir- 
cles stand for the numerical data, and solid curve is the the- 
oretical expression (14). 



IV. WIGNER BAND RANDOM MATRICES 

It is instructive to apply the above analysis to the so- 
called Wigner band random matrices (WBRM). These 
matrices are very useful for understanding generic prop- 
erties of realistic physical systems of interacting particles. 
The WBRM model is described by the Hamiltonian (1) 
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which consists of two parts, one of which is a diagonal 
matrix Hq with increasing entries ej. This part can be 
treated as the "mean field" part of the total Hamiltonian 
H. Another part is a banded matrix Vtj, which is as- 
sociated with the interaction between unperturbed basis 
states. Thus, the model has the following form, 

Hij = ejd,,+V,,, (18) 

where 6ij is the delta-function. It is assumed that ran- 
dom values Cj with the mean spacing D are reordered in 
an increasing way, ej+i > £j. As for the off-diagonal ma- 
trix elements Vij, they are distributed according to the 
gaussian distribution (with the zero mean, < Vij >= 0, 
and the variance < >= Vq) for the matrix elements 
inside the finite band \i — j\ < 6/2, and zero otherwise. 

These matrices have been introduced by Wigner in 
Ref. [20] in application to nuclear physics. A particu- 
lar interest was the form of the strength function in de- 
pendence on the strength of the interaction V. It was 
shown that the form of the SF is the BW for a moderate 
(non-perturbative) strength V, and the semicircle for a 
very strong perturbation. Full analytical treatment of the 
form of the SF for Eq.(18) is given in Ref. [21] with the 
use of the modern approach. In particular, it was found 
that in the transition from the BW to the semicircle, the 
form of the SF is very close to the Gaussian. 

The condition for the SF to be of the Breit- Wigner 
form in the WBRM model can be written as follows [21], 

< Fo < Av; Av = bD = bpo\ (19) 

where the half-width Fq is given by the Fermi golden rule, 

Fo = 27rpon'> (20) 

and Ay is the energy width of the interaction V. 

The left part of the inequalities in Eq.(19) indicates 
the non-perturbative situation for which many of unper- 
turbed basis states are strongly coupled by the interac- 
tion. On the other hand, the interaction should not be 
very strong, namely, the width Fq determined by Eq.( 
20), has to be less than the width Ay of the interac- 
tion in the energy representation. The latter condition is 
generic for systems with finite range of the interaction V. 
One should stress that, strictly speaking, the BW form 
of the strength function is not correct in physical appli- 
cations since its second moment diverges (which assumes 
an infinite range of the interaction). For the WBRM 
model with finite values of b it was shown [20,1] that far 
off the energy range Ay the SF decreases faster than a 
pure exponent. 

In contrast with the TBRI model, in the WBRM model 
the energy scale Ay is well defined that simplifies our fur- 
ther analysis. One can see that instead of the control pa- 
rameter Ay, one can equivalently use the variance A|, of 
the SF, which can be expressed through the off-diagonal 



matrix elements of the interaction, A|; = V^^j for 
i j, therefore, A|; = bV^. As a result, we have 
Ay = 27rA|;/ro and the relation (19) can be written 
in the form, 

£> « Fo < AeV2^. (21) 

Numerical data for the WBRM model have confirmed 
that for Fq ~ 2Ae the form of the SF is quite close to 
the Gaussian. Moreover, it was found that the transition 
from the BW dependence to the Gaussian-like is very 
sharp. Note that the extreme limit of a very strong in- 
teraction, Fq » 2 As, seems to be non-physical, giving 
rise to the semicircle form of the SF. The same effect oc- 
curs for the TBRI model for which the relation V ^ Hq 
means that the residual interaction is much stronger than 
the mean field part. 

Applying the arguments given for the TBRI model, the 
Gaussian decrease of the SF has the form, 

InWoit) = -Alt' = -Fo^ = (22) 

with ^ = t/tc and tc = 2'K/bD. 
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FIG. 4. Return probability Wo{t) for the WBRM model 
(18). Full and open symbols stand for numerical data and 
theoretical expression (22), respectively. 



Numerical data in Fig. 3 obtained for A'' = 924 and dif- 
ferent model parameters, are in a good correspondence 
with our analysis. First, the time dependence of Wo{t) 
on the time scale t is, indeed, of the Gaussian form 
(22). Second, the critical time which divides the two 
characteristic dependencies, nicely corresponds to the ex- 
pression tc ~ 2n/bD « 1.0. As one can see, the transi- 
tion from the Gaussian to the exponential decrease is a 
quite generic property of the wave packet dynamics in 
the systems which are described by the TBRI or WBRM 
models. 
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V. DYNAMICAL MODEL OF A QUANTUM 
COMPUTATION 

Now wc apply our analysis to a physical model of quan- 
tum computation which has no random parameters. Al- 
though a direct application of the results obtained for 
random matrix models to the dynamical models is not 
justified, it is of interest to see whether such a compari- 
son is possible. The model we consider here was recently 
proposed [22] as a simple realization of a solid-state quan- 
tum computation. It describes a one-dimensional chain 
of L interacting 1/2-spins (qubits) that are subject to an 
external magnetic field. In order to have selective reso- 
nant excitation, the time independent part = B^{x) 
of a magnetic field is assumed to have a constant gradient 
along the .x-dircction. This provides different Larmor fre- 
quencies for different spins, = 7„i?^ = + ak, where 
7„ is the spin gyromagnetic ratio and a is proportional 
to the gradient of the constant part of the magnetic field 
(see details in [10]). 

For a specific pulse of the time-dependent part of the 
magnetic field, resulting in a single cubit operation, one 
can derive the time-independent Hamiltonian (1) with 



L-1 



L-1 



fe=0 



k=0 



(23) 



where £,k = n.k [10]. Here the frequency is the Rabi fre- 
quency of the p-th pulse, = (l/2)o"'J'^'^ with o"^'^'^ 
as the Pauli matrices, and = I'^ ± z/^. It is also as- 
sumed that the interaction J between nearest qubits does 
not depend on the indexes k and fc + 1 . 

The unperturbed basis (in which Hq is diagonal) is re- 
ordered according to an increase of the index s which is 
written in the binary representation, s = il-i, iL-2, io 
(with is = or 1, depending on whether the single- 
particle state of the i— th qubit is the ground state or the 
excited one). Therefore, the parameter corresponds 
to a non-diagonal coupling, thus, determining the matrix 
elements Vkn = Vnk = with n ^ k. As one can 

see, in contrast with the TBRI model discussed above, in 
the chosen representation the interaction between parti- 
cles is absorbed by i?o, and V describes the coupling to 
the external magnetic field. 

The problem studied in Refs. [10] was the analysis of 
whether the inter-qubit interaction, as well as the inter- 
action of qubits with the external magnetic field, can be a 
source of a kind of internal decoherence caused by quan- 
tum chaos. It is a wide-spread concern that for many 
interacting qubits the onset of quantum chaos may oc- 
cur even for a very weak interaction, see, for instance, 
[6,23]. This expectation is based on the fact that with an 
increase of number of qubits the level density of many- 
body states increases drastically, thus strongly enhanc- 



ing the delocalization effects due to the interaction be- 
tween qubits. The simple estimate [6] shows that gener- 
ically the threshold for the onset of chaos decreases as 
Jcr ^ 1/L where is a critical inter-qubit interaction 
above which the eigenstates are extended over unper- 
turbed many-body states. 

However, as was found in Refs. [10] , in the case of an 
external magnetic field with a constant gradient along the 
chain of qubits, the onset of chaos is strongly suppressed. 
Although in this case the total energy width may be large 
(proportional to L), the estimate [24] shows a feasibility 
of an experimental realization for a quite large number of 
qubits. Another principally different scheme that allows 
to avoid strong delocalization/chaos effects, is suggested 
and analyzed in Refs. [11] for a quantum computer based 
on electrons on helium. 

Our interest below is to see whether the results for 
the return probability, obtained above in terms of ran- 
dom matrix models, can be applied to the dynamical 
model (23). Note that a strong decrease of the return 
probability Wo{t) can serve as an effective measure of a 
stability of a quantum computation. It can be shown 
that the delocalization effects are directly governed by 
the time-dependence of Wo{t). Strictly speaking, the 
effectiveness of a quantum computation should be an- 
alyzed for the time-dependent model with a large num- 
ber of pulses. However, if for a single pulse the effects 
of quantum chaos are strong, they will be generically en- 
hanced in the presence of many pulses. This is why below 
we restrict ourselves by a consideration of the simplified 
time-independent Hamiltonian (23). We hope that our 
analysis may have also an interest in view of general prob- 
lems of the dynamics of systems with a large number of 
interacting spins. 

First, we start with the so-called non-selective regime 
which is defined by the conditions, CI ^ a ^ J. This 
inequality provides the simplest way to prepare a homo- 
geneous superposition of 2^ many-body states, which is 
needed to start with the implementation of the Slior's or 
Grover's algorithm. The analytical and numerical treat- 
ment of the model (23) in this regime has revealed [10] 
that the constant gradient magnetic field (with a 7^ 0) 
strongly reduces unwanted effects of quantum chaos. 
Specifically, it was shown that in this case the chaos 
border turns out to be independent of the number L of 
qubits, in contrast to the models thoroughly studied be- 
fore [23]. In particular, the quantum chaos may occur 
only for a large coupling and strong interaction J be- 
tween qubits. Another new effect which was found in 
Refs. [10], is that the border of quantum chaos does not 
coincide with the border of delocalization. This pecu- 
liarity is important in view of applications to integrable 
or nearly integrable models for which the quantum chaos 
is absent, however, the delocalization effects can be very 
strong. 
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As was shown in Rcfs. [10] in the non-sclcctive regime 
with n = 100, the delocahzation effects start to play es- 
sential role for J > J,, k, 10. On the other hand, the 
quantum chaos effects, siich as the Wigncr-Dyson distri- 
bution for the spectra statistics, occurs for J > Jg 100. 
Since in this case both the inter-qubit interaction and the 
perturbation due to the magnetic field arc strong, one can 
analyze the dynamics of the return probability Wq (t) in 
connection with the previously discussed results. 
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FIG. 5. Return probability Wo{t) for the dynamical model 
(23) of a quantum computation for L = 8, fi = 100, a = 1 and 
different values of J. Solid curve is the theoretical expression 
(22) for the Gaussian dependence, and full squares, triangles, 
circles are numerical data. 



Numerical data presented in Fig. 5 show a strong de- 
crease of Wo(i) for three values of J. As one can see, 
in all cases the return probability decreases similar to 
the Gaussian form of Eq.(14). The correspondence with 
theoretical predictions should be treated as a good one, 
taking into account the dynamical character of the con- 
sidered model. It is important that the perturbation CI 
due to the magnetic field is strong, therefore, the results 
are almost insensitive to the inter-qubit interaction in a 
very wide region of J. 

Now wc analyze the regime of selective excitation which 
is characterized by the following range of parameters, 
<C J <C a. In this regime each pulse acts selectively on 
a chosen qubit thus resulting in a resonant transition. In 
Ref. [25] this regime was analyzed in detail with the main 
interest to the fidelity of some quantum protocol (for the 
time-dependent Hamiltonian with many pulses). It was 
shown that in this regime for a relatively large gradient 
of the magnetic field there is no any danger of quantum 
chaos, and the dynamics can be analyzed on the base of 
the perturbation theory. It is instructive now to see how 
the return probability Wo{t) corresponds to these results. 
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0.2- 
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FIG. 6. Return probability Wo at time ts ~ 3.0 for the 
model of quantum computation with L = 8, = 0.118, J = 1, 
in dependence on the parameter a which is proportional to the 
gradient of the magnetic field. 

In Fig. 6 we present numerical data demonstrating the 
dependence of Wq for some time ts ~ 3.0. First, one 
can see that for small values of a the perturbation turns 
out to be very strong, leading to a strong decrease of 
Wq. This result demonstrates that even a small inter- 
action between qubits gives rise to a strong leakage of 
the probability from an initially excited state. The ori- 
gin of this phenomenon is the degeneracy (for a = 0) or 
quasi-degeneracy (for small o). With an increase of the 
gradient of the magnetic field, the dynamics turns out to 
be very stable, as is manifestcid by the values of Wo(ts) 
close to one. The analytical estimate of the critical value 
of the inter-qubit coupling above which the fidelity is 
very high, was found [25] to be Ucr ~ 4J. Our data, 
however, clearly show a slightly different value Ucr ~ J- 
A more careful inspection of the numerical data of Ref. 
[25] confirms that, indeed, the fidelity starts to decrease 
at smaller than a ~ 4 values. Our data in Fig. 5 also 
reveal a quite interesting resonance effect that occurs for 
specific values a = 1 and a = 2. This effect has the same 
origin as that found in Ref. [25]. 

In order to see more clearly the role of the gradi- 
ent magnetic field, we performed an additional check of 
the time dependence of the return probability for small, 
a = 0.05, and large, a = 3.0, values of a. The results 
in Fig. 7 manifest a very different behavior of the return 
probability for these two representative values. If for the 
small value of a the decrease of Wo{t) is very fast and 
seems to be non-rccurrcnt, for large a = 3.0 the return 
probability remains very close to one and shows a clear 
recurrence. It is interesting to note that for a = 0.05 
one can see two regimes, the Gaussian one (for t < 2.0), 
and the exponential one (for t > 3.0). Note that due 
to a specific character of the selective excitation (quasi- 
degeneracy for o = 0.05 and very weak perturbation for 
a = 3.0), the comparison with the discussed above ana- 
lytical estimates in this case is not valid. 
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FIG. 7. Return probability Wa{t) for small, a = 0.05, and 
large, a = 3.0, values of the parameter a, for the parameters 
of Fig.6. 



VI. CONCLUSION 



is purely dynamical one (without any random parame- 
ters) . As is now well understood, the effects of quantum 
chaos and/or delocahzation may also arise in the models 
of quantum computation, and lead to a kind of internal 
decoherence. Our numerical data show that in the region 
of parameters that corresponds to a strong delocahzation 
and quantum chaos, the return probability has the Gaus- 
sian form of the decay for a quite a long time, similar to 
what occurs in the random matrix models. Additional 
study of the role of a non-zero gradient magnetic field 
has shown that the fast decrease of the return probabil- 
ity corresponds to earlier results for the fidelity of some 
quantum protocol with many pulses. Therefore, the re- 
turn probability can be considered as a quite sensitive 
quantity for establishing the region of parameters where 
one can expect stable opcrability of quantum computers. 

This research was partially supported by the CONA- 
CYT (Mexico) grant No 43730. 



In conclusion, we have studied the return probability 
Wo{t) in throe models of strongly interacting particles. 
The first model is the model with random two-body inter- 
action, typically used to describe the many-body systems 
of Fermi-particles such as heavy nuclei, many-electron 
atoms, quantum dots, etc. As usual, the two-body ran- 
dom matrix elements are assumed to be gaussian random 
entries, which is known to be reasonable when the inter- 
action is strong and has a complicated form. The main 
result of our analytical and numerical study is that with 
an increase of the interaction, the Gaussian decrease of 
the return probability can last for a long time. The typ- 
ical picture is the following: on some time scale t < tc 
the decrease of Wo{t) is the Gaussian, and for t > tc 'it 
is the exponential one. The Gaussian decrease on the 
scale Q <t <tc can be either weak or strong, depending 
whether the interaction is strong or very strong. Thus, 
the standard exponential decrease, associated with the 
Fermi golden rule, is not correct for a very strong per- 
turbation. Numerical data confirm the analytical predic- 
tions. 

In order to elucidate the meaning of our analytical es- 
timates, we have also analyzed the model of Wigner band 
random matrices which captures essential features of the 
chaotic systems of interacting particles. This model turns 
out to be very effective since it demonstrates in a very 
transparent way the dependence of the dynamics on few 
global parameters of physical significance. Our numer- 
ical data reflect the generic properties of the interplay 
between the Gaussian and exponential decrease of the 
return probability. 

Finally, we have analyzed the model which was recently 
proposed as an implementation of a quantum computer. 
In contrast with the random matrix models, this model 
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